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Abstract. Differential equations arising in fluid mechanics are usually derived from 
the intrinsic properties of mechanical systems, in the form of conservation laws, and 
bear symmetries, which are not generally preserved by a finite difference approxi- 
mation, and leading to inaccurate numerical results. This paper develops a method 
that enables us to build a scheme that preserves those symmetries. The method is 
based on the concept of the differential approximation. A comparison of numerical 
performance of the invariant schemes, standard ones and higher order one has been 
realised for the Burgers equation. 

1. Introduction. Lie groups were introduced by Sophus Lie in 1870 in order to 
study the symmetries of differential equations, yielding thus analytical solutions. 
Literature provides substantial works and applications, [3J, [!]• Symmetry groups 
can be determined by an automatic procedure, but it often turn out to be tedious 
and induce errors. A large amount of packages using symbolic manipulations of 
mathematical expressions have been written. We mention here some of those works: 
Schwartz [T7], Vu and Carminati|14j. Herod [T5], Baumann [If)], Cantwell [5]. 
In this paper we are interesting in the application of the theory of Lie group to 
numerical analysis. 

Finite difference equations used to approximate the solutions of a differential equa- 
tion generally do not respect the symmetries of the original equation, and can lead to 
inaccurate numerical results. Various techniques, that enable us to build a scheme 
preserving the symmetries of the original differential equation, have been studied. 
One of these techniques consists in constructing an invariant scheme from a given 
one by applying the method of the moving frame in [7] , [8] . Another one consists in 
constructing an invariant scheme with the help of the discret invariants of its sym- 
metry group [9] , [TU] , [IT] , [T2] , [13] and provides the building of symmetry-adapted 
meshes, in preserving the differential equation symmetries. This technique is based 
on a direct study of the symmetries of difference equations and lattices. 
Yanenko [2] and Shokin [1 , proposed to apply the Lie group theory to finite dif- 
ference equations by means of the differential approximation. Thus, they have set 
down conditions under which the differential representation of a finite difference 
scheme preserves the group of continuous symmetries of the original differential 
equation. They provide a dissipative scheme, which is called invariant scheme. 
The resulting scheme is independent of any change of the reference frame, and its 
differential representation is invariant under the symmetries of the original equa- 
tion. Ames, Postell and Adams [6j have already used the approach of Yanenko 
and Shokin to present invariant schemes in which terms are added to the original 
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difference scheme. They showed that, in specific cases, the invariant scheme is as 
accurate as high order numerical methods. 

In this paper, we focus on the last approach. The method is implemented on some 
standard schemes for the Burgers equation. A comparison is made between the 
numerical solutions of these schemes and the invariant scheme. 
The paper is organized as follows. Definitions and invariance condition for differen- 
tial equations are provided in section 2. Section 3 recalls the approach of Yanenko 
and Shokin. Section 4 concentrates on classical schemes. In section 5, we present a 
method that enables us to build an invariant scheme with respect to an otherwise 
lost symmetry. 



2. Definitions and invariance condition for differential equations. A r- 

parameter Lie group G r of point transformations in the Euclidean space £{x,u) 
can be written under the form: 

G r = {x* = (pi(x,u,a); u* = ipj(x,u,a), i = 1, . . . , m; j = 1, , n} (1) 

Consider a system of l th -oidei differential equations: 

F x (x,u,u ( - kl \u<- kl ' k2 \...,u( kl ~ k ri) = 0, X = l,...,q (2) 

Denote by u i k ^--- k p) the vector, the components of which are partial derivatives of 
order p, namely, u f 1 '" kp) = dx ^^ 3 dXk ^ j = l,...,n and kx,...,k p S {l,...,m}. 
Denote by x = (x\, . . . , x m ) the independent variables, u = (ui, . . . , u n ) the depen- 
dent variables, and [xk t ■ ■ ■ Xk p ) a set of elements of the independent variables. 
Equation (J2]) is a subset of the Euclidean space £Cx,u, u^ kl \ . . . , y^ kl '" kl '). In order 
to take into account the derivative terms involved in the differential equation, the 
action of the group G r of transformations in the space £{x, u) needs to be extended 
to the space of the derivatives of the dependent variables. 

Denote by G r a r-parameter Lie group of point transformation in the space 
£ (x, u, v,( kl \ . . . , u( fc i— fc 0) of the independent variables, dependent variables and 
the derivative of the dependent variables with respect to the independent ones. 
The /^-prolongation operator of G r is: 

i = 1, . . . , m; j = 1, . . . ,n; a = 1, . . . , r. 

, T]" , a - A ' and ' are given by: 

1 da a \ a=o 3 da a \ a—o 3 T^x kl & x i T^ x ki 

^■■■M = ^ £ ^L, = 2 , 

3 ' Dx k ~[dxtdx kl . . .dx ko l Vx ko 

n V d v-^ dv,j d 

where: = h > — 

T>x k dx k |— J dx k duj 

The system of l th -order differential equations is invariant under the group G r if 



and only if: 



-L^r x \ =0, a = l,...,r; A = l,...,o (4) 



LIE GROUP STABILITY OF FINITE DIFFERENCE SCHEMES 



3 



3. Lie group for the differential approximation. The finite difference scheme, 
which approximates the differential system l[2|). can be written as: 

A x (x,u, h,Tu) = 0, X = l,...,q (5) 

where h = (hi, ft,2, . . . , h m ) denotes the space step vector, and T — (T±, T2, . . . , T m ) 
the shift-operator along the axis of the independent variables, defined by: 

T i [u](x 1 ,x 2 , ■ ■ ■ , Xi-i,Xi,Xi+i,. . ■ ,x m ) = u(xi,x 2 , . . . , Xi-i, Xi + hi,x i+1 , . . . ,x m ). (6) 
Definition 1. The differential equation: 

V x (x,u,u^ 1 >,..., n( fe i- fe i')) = T x (x,u,u^\...,u^- k ^) 

s m 

+ E X>i)^(w (fel) ,. ■ ..« ( 1 "^J, 
0=1 i=l 

A = 1, . . . , q; I' = max (x>i) l' x ti (7) 

is called the s'^-order differential approximation of the finite difference scheme 
([5]). In the specific case s — 1, the above equation is called the first differential 



approximation. 



Equation ([7]) is obtained from equation ^ by applying Taylor series expansion 
to the components of Tu about the point x — (x%, . . . , x m ) and truncating the 
expansion to a given finite order. Denote by G' r a group of transformations in the 
space £(x, u, h): 

G' r = {x* = (f>i (x, u, a); u* = ipj (x, u, a); h* = ipi (x, u, h, a), i = 1, . . . , m; j = 1, . . . , n} (8) 
by L Q ' the basis infinitesimal operator of G' r : 

L a ' = €f(x,*)^-+i7f(x,«)^-+<f(x I u,h)^- ) a = l,...,r (9) 
oxi ouj ahi 

where 

Cf = |^| , a = l,...,r (10) 
aa a la— 

and by G'r ' a group of transformation in the space £ h,u( kl \ . . . , u^ 1 *"*'')). 
The ^^-prolongation operator of G* r , La ^ can be written as: 



71 I Q 

l(0 = w + EE<t ;^-M (ID 

j=ip=i 9«j 

Theorem 1. T/ie differential approximation ^ is invariant under the group Gr 
if and only if 

h^" > V X ((x,u,u ( - kl '>,...,u (kl - k i' ) )\ x Q =0, a=l,...,r; \ = l,...,q (12) 

or 

s m 

[Eg ) J rA +£g ) (X)$3(Ai)'''^)]L =0, a = l,...,r; A = l,..., 5 (13) 
0=1 «=1 73 ° 

Theorem [T] provides the equations which enable us to obtain the symmetry groups 
of the differential approximation. The unknowns are the infinitesimal functions 77" , 
£f and i = 1, . . . , m; j = 1, . . . , n. The infinitesimals a ^ kl '— > feo ) j j = l ; . . . ; ri, 
are functions of the partial derivatives of 77" and £° . 

Equation (fT3")l is simplified by means of the condition ([7]) . They lead to an overdeter- 
mined system of differential equations, the unknowns of which are the infinitesimal 
functions. 
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4. The specific case of the Burgers equation. 

4.1. Symmetries of the Burgers equation. The Burgers equation can be writ- 
ten as: 

F{x,t,u, v, u x ,ut,Uxx) = u t + u u x - v u xx = (14) 
where v > is the dynamic viscosity. 

Denote by G a group of transformations of the Burgers equation in the space 
£(x,t,u, v) of the independent variables the dependent variable u, and the 

viscosity v. The viscosity is taken as a symmetry variable in order to enable us to 
take into account variations of the Reynolds number. 
G is a set of transformations acting smoothly on the space £(x, t, u, v). 
The six-dimensional Lie algebra of the group G is generated by the following oper- 
ators: 

d d a a a 

Li = — , L.2 = — , L3 = x h 2t u — 

dx dt dx dt du 

9 9 a . a a a a a a , . 

L 4 = xt ht h(-ut + x) — , L 5 = t 1 , he = -t \-u Yu — (15) 

dx dt du dx du dt du du 

which respectively correspond to: 

• the space translation : (x,t, u, v) 1 — > (x + ei,t, u, v); 

• the time translation : (x,t,u, u) 1 — ► (x,t + €2,u, v)\ 

• the dilatation : (x,t,u,u) 1 — ► (533;, e|t, £3 1 u, v); 

• the projective transformation : (x,t,u,u) 1 — > f 1 _^ 4t , 1 _\ it , ^<?4 +u(l - e<ii),v\; 

• the Galilean transformation : {x,t,u, v) 1 — ► (x + e 5 t,t,u + es,u); 

• the dilatation : (x,t,u,u) 1 — > (x, eg 1 t, eeu, eeu). 

(ei)i=i g are constants. 

4.2. Symmetries of first differential approximations. Denote by h the mesh 
size, r the time step, N x the number of mesh points, N t the number of time steps, 
and i € {0, . . . , N t }, n 6 {0, . . . , N x } the discrete approximation of u(ih, nr). 
In order to shorten the size of the finite difference scheme expressions, we use the 
following notations introduced by Hildebrand in [18]: 

U" ! - U n _ 1 U n 1 + U n _ 1 



<> K) = 1 1 K) = 

h h 

The Burgers equation can be discretized by means of: 

• the FTCS (forward-time and centered-space) scheme: 

< +1 -< + ^r-4<=o 

h v 9 '» h2 1 



• the Lax-Wendroff scheme: 



4< 



where: 



2h 2 L 1 v 2 y » v 2 ' l J 2 L/1 4 * 



)/)3 L £J«? 5 2 (£i<) - + v -^[^{\) 

• the Crank-Nicolson scheme: 
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Linear stability properties and the related orders of approximation are displayed in 
Table [T] (where CFL = <f, S = % and S* = {v + 9 ^P ± )jp)- 



Scheme 


Stablility condition 


Error 


FTCS 


S<i, CFL < 1 


0(t, h 2 ) 


Lax-Wendroff 


S* < i, CFL < 1 


0(r 2 ,h 2 ) 


Crank-Nicolson 


unconditional stability 


0(r 2 ,h 2 ) 



Table 1 . Table of finite difference schemes 



Consider m n as a function of the time step r, and of the mesh size h, expand it at 
a given order by means of its Taylor series, and neglect the o(t q ) and o(h^) terms, 
where a and (3 depend on the order of the schemes. This yields the differential 
representation of the finite difference equation. 
The following differential representations are obtained: 

• for the FTCS scheme: 

1 2 n r ^ 2 2\ t/ ^ 2 

Uf + -(« )x — V Uxx + -Q2 + —{U ) x xx — U xxxx = 

£t L i£i 

• for the Lax-Wendroff scheme: 

l , 2, T 2 h2 , 2^ 

Ut + -(« )x — V U xx + —93 + TTTV" )xxx — — Uxxxx = 

^ O L iu _L _ 

• for the Crank-Nicolson scheme: 

1 n 2/93 1/2 v / \ \ h 2 . 2\ u h 2 
Ut + -(u )x - vuxx + t ( — + -{g 1 + ug 2 )x - -(g2)xx) + TTTl" )xxx — Uxxxx = 

^ O 4t rr L^ L^ 

where S i = -(^ + uu xx , 92 = ( - 91^ + v(gi) xx , 93 = ( - g 2 u - g 2 ) x + * (g 2 ) xx 

Denote by G the group of transformations of a first differential approximation in 
the space E{x,t,u,h,T,v) of the independent variables and the dependent 

variable u, the step size variables (h, r) and the viscosity v. 
The ^'^-prolongation of G' can be written as: 

t£"> = gf— 1 - 1 t;" ° I y^ CT ".( fc i- fc p) 8 g 8 , g 9 , r 8 (16) 

9ie cM <9m J du^ kl "' kp ^ &h & T ® v 

where I' has been defined in definition [TJ 

Theorem[TJenables us to obtain the necessary and sufficient condition of invariance 
of the first differential approximation V: 

= (17) 

\v=o 

Theorem Q] is applied to the differential representations of the above schemes. 
The resolution of the determining equations of each first differential approximation 
yields the 4-parameter group: 

gf=a + 6x, ff = c+(2b-d) t, r) a = (-b + d)u (18) 

£f = 6 h, CS = ( 2b - d ) T : ° a = eV 

The 4-dimensional Lie algebra of G 1 is generated by: 

a T a ,, a a a a a 

Li = — , L 2 = — , L' 3 =3! h 2* u h/i \-2t — 

dx dt dx dt du dh dr 

• a a a a , , 

L' 4 = -t hw t Yv — (19) 

dt du dr du 

These operators are respectively related to: 
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• the space translation : (x, t, u, h, t, v) i — ► (x + ei,t,u, h,T,u); 

• the time translation : (x,t, u, h,r, v) i — ► (x,t + e2,u, h,r, v)\ 

• the dilatation : (x, t, u, h, r, u) i — > (e 3 x, e§t, eg 1 u, £3/1, e|r, v)\ 

• the dilatation : (x, t, U, h, T, u) i > (x,eJ 1 t,e4U,h,eJ 1 T,e i u)] 

where (ei)i=i,...,4 are constants. 

The above finite difference equations are preserved by the space translation, the 
time translation and both dilatations. 

Approximating the Burgers equation by the above finite difference equations results 
in the loss of the projective and Galilean transformations. 

5. The invariant scheme. 

5.1. Invariant scheme construction. An invariant scheme is constructed in such 
a way that the related differential approximation preserves the symmetries of the 
Burgers equation. We propose to approximate the Burgers equation by the following 
finite difference scheme: 

< +1 - + i / s _ MS A» _ * rp - £t) «) - f n" w g+ - n? yS-)tff = o (20) 

t h y 6 A 2'" h 2 v 12 M w V »+5 »-5 J 1 v y 

where f2™ = t n , u^) is defined next so that the related differential represen- 
tation is preserved by the symmetries of the Burgers equation. The scheme has 
second-order accuracy in space and first-order accuracy in time. The derivatives 

\\) x and u xx are approximated by fourth order accuracy difference expressions: 

(f-^)(<) = (u x -^u 5 ^ + o( h % { ^-^ ){<) = {uxx -^ U(>x y; + o( h «) (21) 

The truncation error of the difference scheme (|20j) can be written as: 

6 = ^-u tt -h 2 (nu x ) +0(T 2 ) + 0(h 4 ) 
2 V / x 

u u is replaced by an expression involving partial derivatives with respect to x, by 
using the Burgers equation: 

u 2 

Uu = {u 2 u x ) x - v{uu xx ) x - <j(—) xxx + v 2 u xxxx (22) 

Replacing the previous expression in the truncation error leads to: 



where C = —u 2 — h 2 il. 
2 

It is convenient for the calculation of C that the truncation error is reduced to: 

6 = (Cu x ^+0{T 2 ) + 0(h 4 ) 

The related finite difference scheme is the following first order accuracy in time 
and second order accuracy in space: 

™ +1 -- n 1 fiS 3 .,u 2 . n 1 , r2 <5 4 



us - ^-) (— )" - uA-is 2 - —)(«?) - (n» x Ai - n n 1 a_iW 



t /i v 6 ' v 2 '« /i 2 v 12 ' 

vt ( _ M<5 2 , n . „ A«5 2 ,„ A v 2 T 5 A vr ^ ,u 2 . n 

H M , 1 rT" ( u - , 1 ) — u - 1 ;T \ u ■ 1 ) 7 u i H 5" ( ) . = U (2o) 

2 V *+± ft 2 v t-J h 2 ^ '-iV 2 h 4 1 2 V 2 ; * V ; 

and the differential approximation can be written as: 

P(:r, t, «, Uo;,u«,Uaia!) = u* + u u x — v u xx + (Cu x ) x = (24) 

The von Neumann stability analysis of scheme (|23p under a linearized form provides 
the following necessary conditions for S, CFL and fl T — fir: 

45 1 

CFL 2 - 2S - 2Q T < 0, < 2S 2 + Q T < - (25) 

~ 3 ~ 2 
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If £1 takes is sufficiently close to zero, these conditions become then sufficient for 
the linear formulation. 

5.2. Calculation of the artificial viscosity term. Here we describe the method 
for determining the artificial viscosity term (Cu x ) Xl which is constructed in such 
a way that the differential approximation (|24p is preserved by the symmetries 
of the Burgers equation. C is a function of the variables (x,t,u,T,h), and also 
depends on the partial derivatives of u with respect to x: u x and u xx . C = 
C(x,t,h,T,u,u x ,u X x)- The necessary and sufficient condition for the differential 
approximation to be an invariant of the Burgers equation symmetry group is: 

{u t + u u x - u u xx )\ + l i <? ) ((Cu x ) x )\ =0 (26) 

\v=o \v=o 
Equation (|26|) provides the determining equations of the symmetry group of equa- 
tion (|24p . The determining equations involve partial derivatives of the unknown 
function C and partial derivatives of the infinitesimal functions of G , which is the 
symmetry group of the differential representation of the invariant scheme. 
The infinitesimal functions of G have the following expressions: 

£f = a + b x + c t + d tx, ($ = e + d t 2 + (2b - f) t, Cf = b h, (27) 

Ca=(2b-/)r, v a = c + d x + (-b - d t + f) u, 8 a = fv 

The determining equation with respect to the unknown function C is simplified in 
using the infinitesimal functions of each subgroup of G . 

The determining equations of each subgroup of G provides the following linear 
partial differential equations and the expressions for C: 

• the space translation -^-C = =>• C = Ci(t, h,r,u,u x ,u xx ); 

• the time translation J^C = => C = Cz{x, h,r,u, u x ,u xx ); 

• the dilatation x^C+2t§ l C-u-^C+h-§ K C+2T-§^C = 0=> C = C 3 (4f, |,nx, ^,u x ,u x . 

• the projective transformation J-C = 0, J-C = 0, = 0, t 2 4c + 2-^-C = 



C = C 4 (h, 



• the Galilean transformation -S-C + t-§-C =*• C = Cs( t t,h,T,u,u x ,u xx ); 



du dx 

• the dilatation -t-S-C + u-S-C + u-g-C - rg-C = => C = \C & (x, h, J). 

ot ou ou or t u v ' ' * ' 



5.3. Numerical application. The numerical resolution of the Burgers equation 
has been implemented for scheme (|23p , the standard schemes (cf. section 14.21) 
and a scheme with second-order accuracy in time and fourth-order accuracy in 
space, which is obtained from the invariant scheme when C = 0. The solutions 
are calculated in the reference frame (Fl) and in the one (F 2) resulting from the 
Galilean transformations (x, t, u, v) i — ► (x + t, t, u + 1, v). The artificial viscosity 
has the following expression: 

C = -0.01t(tu - x) 2 (u x ) 2 (28) 

1 r 

SI = J2^~2 u2 ~ @) 1S m a sufficiently small neighborhood of zero that we have the 

sufficiency of conditions (|25|) for the linear formulation. 

The problem consists in solving the following differential system: 

u t + uu x - vu xx = 0, x e [0, 40], t 6 [0, 20] 

«(x,0) = /(x), w(0,t) = s(t), u(40,t) = h(t) 

The initial and boundary conditions, /, h, and g are provided by an exact solution 
of the Burgers equation: 

, \ (x - 2t)/(t + 0.1) , , 

V ' ; 1 + u 2 Vt + 0.1 exp ((x - 2t) 2 /(4^(t + 0.1))) V ' 

Figures [1] [3] and [5] show the time evolution of the L 2 -norm of the error for the 
considered schemes, for specific values of the CFL number and the mesh Reynolds 
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number Ren- Figures [2j |4] and [6] display the variations, as functions of the space 
variable, of the numerical solutions of the considered schemes for the specific value 
t = 5. In each frame, the numerical solutions are compared to the exact one. 
The error analysis of the invariant scheme in the reference frame through the features 
of the truncation error and the graphical representation of the norms of the error 
(cf. Figures [TJ [3] and [5|) allows to say that the invariant scheme is dissipative and 
slightly dispersive. 

The presence of the dissipative term {Cu x ) x in the differential representation of 
the invariant scheme and the presence of the higher order error terms involving 
the even-order derivative uq x (cf. Equation (|21j) ) show that the scheme produces 
numerical damping. Particularly, the amplitudes are not correctly represented for 
high frequencies, since the solution is subjected to rather rough variation during 
the first iterations. The dissipation is stronger for Reh = 2, CFL = 0.08 in the 
reference frame (see Figure [4]). Moreover, the presence of higher order error terms 
involving the odd-order derivative u$ x corresponds to a phase error. 
The non-invariant schemes are more altered by the change of the frame than the 
invariant one. Moreover, the invariant scheme appears to be as accurate as the 
higher order one in the frame (F2). 




FIGURE 1. Evolution of the error L 2 -norm in (Fl) and (F2). Re h = 2, CFL = 0.04 




Figure 2. Space variation of the numerical solutions of the schemes and the 
exact solution in (Fl) and (F2). Re h = 2, CFL = 0.04 
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FIGURE 3. Evolution of the error L 2 -norm in (Fl) and (F2). Re h = 2, CFL = 0.08 




Figure 4. Space variation of the numerical solutions of the schemes and the 
exact solution in (Fl) and {FT). Re h = 2, CFL = 0.08 
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Figure 6. Space variation of the schemes numerical and of the exact solution 
in (Fl) and (F2). Re h = 3, CFL = 0.08 



